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ABSTRACT 

We present extensive radio observations of the nearby Type Ibc supernovae (SNe Ibc) 2004cc, 2004dk, and 
2004gq spanning At « 8- 1900 days after explosion. Using a dynamical model developed for synchrotron 
emission from a slightly decelerated Shockwave, we estimate the velocity and energy of the fastest ejecta and 
the density profile of the circumstellar medium. The Shockwaves of all three supernovae are characterized by 
non-relativistic velocities of v » (0.1 — 25)c and associated energies of E « (2— 10) x 10 47 erg, in line with the 
expectations for a typical homologous explosion. Smooth circumstellar density profiles are indicated by the early 
radio data and we estimate the progenitor mass loss rates to be M sa (0.6- 13) x 10~ 5 M© yr" 1 (wind velocity, 
v K = 10 3 km s" 1 ). These estimates approach the saturation limit (M rs 10 -4 M yr" 1 ) for line-driven winds from 
Wolf-Rayet stars, the favored progenitors of SNe Ibc including those associated with long-duration gamma-ray 
bursts. Intriguingly, at later epochs all three supernovae show evidence for abrupt radio variability that we attribute 
to large density modulations (factor of ~ 3-6) at circumstellar radii of (1-50) x 10 16 cm. If due to variable 
mass loss, these modulations are associated with progenitor activity on a timescale of ~ 10-100 years before 
explosion. We consider these results in the context of variable mass loss mechanisms including wind clumping, 
metallicity-independent continuum-driven ejections, and binary-induced modulations. It may also be possible that 
the SN Shockwaves are dynamically interacting with wind termination shocks, however, this requires the environ- 
ment to be highly pressurized and/or the progenitor to be rapidly rotating prior to explosion. The proximity of the 
density modulations to the explosion sites may suggest a synchronization between unusual progenitor mass loss 
and the SN explosion, reminiscent of Type Iln supernovae. This study underscores the utility of radio observations 
for tracing the final evolutionary stage(s) of SN progenitor systems. 

Subject headings: supernovae 



1. INTRODUCTION 

Type Ibc supernovae (SNe Ibc) are commonly understood to 
mark the death of massive stars, M > 10 M Q . A relatively rare 
sub-group, SNe Ibc represent just 10-20% of all local SN dis- 
coveries (Cappellaro, Evans & Turatto 1999; Smartt 2009; Li et 
al. 201 1). They are distinguished from other core-collapse SNe 
by the lack of hydrogen features in their optical spectra (see Fil- 
ippenko 1997 for a review) indicating that their massive hydro- 
gen envelopes are removed prior to explosion (Eliasef al. 1985; 
Filippenko & Sargent 1985; Wheeler & Levreault 1985; Uo- 
moto & Kirshner 1985; Clocchiatti et al. 1986; Woosley, Heger 
& Weaver 2002). The lack of helium absorption features may 
further divide the sample into SNe Ic (helium-poor) and SNe lb 
(He-rich; Matheson et al. 2001 but see Hamuy et al. 2002). 

The physical mechanism and timescale on which SNe Ibc 
progenitors lose their envelopes, however, remains an open 
question. In a favored progenitor model, the metallicity- 
dependent (and often clumpy) line-driven wind of an iso- 
lated Wolf-Rayet (WR) star expels its own hydrogen enve- 
lope (Begelman & Sarazin 1986; Woosley, Langer & Weaver 
1995). However, as argued by Smith & Owocki (2006), 
metallicity-independent continuum-driven winds and/or hydro- 
dynamic eruptions may play a key role in the formation of 
Wolf-Rayet stars. A binary mechanism such as Roche lobe 
overflow, accretion or a common envelope phase could also re- 
move the envelope of a massive progenitor (Wheeler & Lev- 



reault 1985; Podsiadlowski, Joss & Hsu 1992; Yoon, Woosley 
& Langer 2010). Both single- and binary-star progenitors 
are similarly theorized for the ~ 1% of SNe Ibc associated 
with long-duration gamma-ray bursts (GRBs), characterized by 
a central engine-driven relativistic outflow (e.g., MacFadyen, 
Woosley & Heger 2001; Podsiadlowski, et al. 2004). How- 
ever, the critical difference between ordinary SNe Ibc progen- 
itors and those of GRB-SNe remains debated (e.g. Fryer et al. 
2007). 

Recent studies of the host galaxies of SNe Ibc and GRB- 
SNe have revealed surprising clues. Specifically, explosion site 
metallicity measurements suggest that GRB-SNe favor metal- 
poor environments more often than SNe Ic (Modjaz et al. 2008; 
Levesqueef al. 2010). In parallel efforts, it has also been re- 
ported that SNe Ic reside in higher metallicity environs than 
SNe lb (Anderson et al. 2010; Modjaz et al. 201 1). These ap- 
parent environmental differences have been interpreted as ev- 
idence for a metallicity-dependent mass loss mechanism (e.g., 
Arcavi et al. 2010), such as line-driven winds (Castor, Abbott 
& Klein 1975; Conti 1978). However, direct observational con- 
straints on the local (sub-parsec) circumstellar environments of 
the explosions are required to test this hypothesis. 

The mass loss histories of single- and binary-star progenitors 
each produce unique density distributions in the circumstellar 
environment (e.g., smooth stellar wind profile, dense common 
envelope shell). The dynamical interaction between the ex- 
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panding SN Shockwave and the circumstellar medium (CSM) 
gives rise to non-thermal radio and X-ray emission that, in 
turn, traces the local mass distribution (Chevalier 1998). Ra- 
dio and X-ray observations of young SNe Ibc are, in fact, the 
only way to reveal the mass loss histories of of their progenitor 
stars in the final evolutionary stage. Motivated thus, the sam- 
ple of SNe Ibc with detailed radio studies is growing, primar- 
ily as part of the Very Large Array Intensive Study of Naked 
Supernovae 3 (VISioNS; Berger, Kulkarni & Chevalier 2002; 
Berger et al. 2003a; Soderberg et al. 2005, 2006a,b,c; Soder- 
berg 2007; Soderberg et al. 2008, 2010a,b). 

The VISioNS study has revealed several key findings to date, 
including (i) SNe Ibc radio luminosities span four orders of 
magnitude, L v « 10 25 - 10 29 erg s" 1 Hz" 1 (Soderberg 2007), 
and (ii) nearly half of all radio-detected SNe Ibc show evi- 
dence for small-scale (factor of ~ 2) light-curve modulations 
(e.g., SN2003bg; Soderberg et al. 2006a). Both of these find- 
ings are attributed to a broad diversity in the mass loss histories 
of their progenitor stars (Chevalier & Fransson 2006). In paral- 
lel efforts, detailed radio studies of GRBs have revealed similar 
diversity in their circumstellar environments (e.g., Chevalier, 
Li & Fransson 2004), however, a direct comparison of the two 
samples has yet to be conducted. 

Here we present extensive radio observations for three Type 
Ibc SNe - 2004cc, 2004dk, and 2004gq - all of which show 
large (factor of <~ 10) flux modulations at late epochs. Through 
our detailed modeling of the radio light-curves, we derive the 
properties of the Shockwave (velocity, energy) and those of 
the circumstellar medium (density profile, progenitor mass loss 
rate). In all three SNe, we attribute the observed flux modula- 
tions to density modulations in the local environment. We com- 
pare the radio properties with those of other SNe Ibc and nearby 
central engine-driven explosions, including GRB-SNe 1998bw 
(Kulkarni et al. 1998) and 2006aj (Soderberg et al. 2006b) and 
the relativistic SN2009bb (Soderberg et al. 2010a). In §2 we 
present the observations and in §3 we describe the dynamical 
model used to extract the properties of the Shockwave and lo- 
cal environment and present the results of modeling the early 
radio data. In §4 we compare the models with the late-time 
radio data characterized by flux modulations. We discuss the 
implications of the observed flux variations in the context of 
circumstellar density modulations in §5 and, finally, in light of 
different channels for progenitor mass loss. 

2. RADIO OBSERVATIONS 

We observed SNe 2004cc, 2004dk, and 2004gq with the Very 
Large Array 4 (VLA) shortly after optical discovery as part of 
the VISioNS survey. In each case we detected a radio source 
coincident with the optical position and subsequently initiated 
an intense VLA follow-up campaign to study the temporal and 
spectral evolution of the radio emission. In Table 1 we sum- 
marize the basic properties of these SNe, including the spectro- 
scopic classification, host galaxy distance estimate, and approx- 
imate explosion date. We adopt host galaxy distance estimates 
and integrated apparent magnitudes from NED 5 , including cos- 
mology independent distances, when available. 

3 Here we use "naked" to describe supernovae from progenitor stars that have 
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Radio data were collected at multiple frequencies between 
1 .43 and 43.3 GHz (Tables 2-4). The data were obtained in var- 
ious array configurations as denoted in the Tables. All VLA 
observations were taken in the standard continuum observing 
mode with a bandwidth of 2 x 50 MHz. At observing fre- 
quencies above 22 GHz, we included reference pointing scans 
to correct for the systematic 10-20 asec pointing errors of the 
VLA antennas. We used primary calibrators 3C48, 3C147, and 
3C286 for flux calibration, while phase referencing was per- 
formed using a calibrator within 10 degrees. Data were reduced 
using standard packages within the Astronomical Image Pro- 
cessing System (AIPS). We fit a Gaussian model to the radio 
SN in each observation to measure the integrated flux density. 

The radio light-curves for the three SNe are shown in Fig- 
ures 1-5, spanning At w 8- 1900 days after explosion, and 
are discussed in detail below. Their peak spectral luminosi- 
ties span L v « (2-6) x 10 27 erg s" 1 Hz" 1 , similar to other SNe 
Ibc and a factor of ~ 100 below the radio afterglow luminosi- 
ties of nearby GRB-SNe and the engine-driven SN 2009bb. In- 
triguingly, all three SNe reveal evidence for large-scale (factor 
of <~ 10) light-curve variations at late time. We estimate the 
"abruptness" of the flux variations by determining the time in- 
terval between the epoch at which the modulation is first seen 
and the epoch of the previous observation (St) and normalize 
by the time since explosion, Af . 

2.1. SN2004cc 

SN 2004cc was discovered on 2004 Jun 10.26 UT by the Lick 
Observatory Supernova Search (LOSS; Li et al. 2000) using 
the Katzmann Automatic Imaging Telescope (KAIT; Monard 
& Li 2004). The host galaxy, NGC 4568, is classified as SAbc 
with an absolute B-band magnitude of M B « -19.6 mag as- 
suming a distance of d s=s 18 Mpc. Spectral analysis revealed 
SN2004cc to be a highly reddened SN Ic (Matheson et al. 
2004; Foleyef al. 2004). An unfiltered light-curve indicates a 
peak optical magnitude of M fa -16.2 mag, less luminous than 
typical SNe Ibc (Li et al. 2011). However, this estimate does 
not account for host galaxy extinction which is likely to be sig- 
nificant (see Drout et al. 2011 for a discussion of extinction 
diagnostics for SNe Ibc). 

Our first VLA observation was carried out on 2004 Jun 
17 UT, roughly At w 15 days after explosion. We moni- 
tored the radio emission from SN 2004cc at frequencies, v = 
8.5, 15.0,22.5 and 43.4 GHz through Af « 1700 days (Table 2 
and Figure 2). At 8.5 GHz, the lightcurve reaches a peak of 
F v w 4.5 mJy at Af w 28 days with a spectral luminosity of 
1.7 x 10 27 erg s" 1 Hz" 1 . An abrupt steepening of the observed 
flux was observed at all frequencies at Af « 40 days after ex- 
plosion with a temporal profile at least as steep as AF„ oc f" 3 4 
(v = 8.46 GHz). Subsequently, the SN re-brightened at Af w 
135 days. We estimate the abruptness of these two episodes to 
be (St /At) « 0.27 and (St /At) £ 1.2, respectively. 

2.2. SN2004dk 

SN2004dk was discovered by KAIT on 2004 Aug 1.19 UT 
(Graham & Li 2004). The host galaxy, NGC 6118 (morphol- 
ogy SAcd), has an absolute magnitude of M B « -19.4 at d ~ 23 

their hydrogen envelope prior to explosion. 

: Foundation operated under cooperative agreement by Associated Universities, 
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Mpc. Although SN 2004dk was originally classified as Type 
Ic, it was re-classified as a Type lb following the detection of 
helium features (Filippenko et al. 2004). A late-time spec- 
trum published by Maeda et al. (2008) reveals an ordinary 
hydrogen-poor SN lb spectrum. Detailed optical photometry 
for SN 2004dk were reported by Drout et al. (201 1) indicating 
a broad light-curve with a peak magnitude of My fa -18.2 mag. 

We first observed the SN with the VLA on 2004 Aug 7 
UT, roughly At fa 8 days after explosion. Table 3 summa- 
rizes our observations of SN 2004dk, taken at frequencies of 
v = 4.9, 8.5, 15.0 and 22.5 GHz. Comparable in spectral lu- 
minosity to SN 2004cc, it reaches a peak of of L v = 1 .6 x 10 27 
erg s" 1 Hz" 1 at Af fa 14 days after explosion in the v = 8.5 GHz 
band. Similar to the case of SN 2004cc, radio monitoring re- 
vealed a re-brightening at Af fa 4 years after explosion (Stock- 
dale et al. 2009) which was not evident in the previous obser- 
vation at Af fa 200 days. We therefore estimate the abruptness 
of the modulation to be <; 7.7. 

We note that X-ray emission was also detected from 
SN2004dk with XMM on 2004 Aug 12 UT (Af fa 13 days) 
with a flux, Fy = 2.7 ± 0.8 x 10" 14 erg cm" 2 s" 1 (0.5-8 keV; Poo- 
ley 2007; Stockdale et al. 2009). At the distance of NGC 61 18, 
this implies an X-ray luminosity of Lx fa 2 x 10 39 erg s" 1 and 
comparable to the X-ray luminosities of other ordinary SNe lb 
observed on a similar timescale (e.g., SN2008D; Soderberg et 
al. 2008). At this epoch, the spectral index (j3; F v oc v&) ob- 
served between the radio and X-ray bands is (3 RX ~ -0.77 ± 
0.03. 

2.3. SN2004gq 

SN 2004gq was discovered independently by KAIT and the 
Stazione Astronomica di Sozzago Supernova Search on 2004 
Dec 1 1.36 UT (Pugh et al. 2004). The host galaxy, NGC 1832 
(morphology SBbc) and has an absolute B-band magnitude of 
M B fa -20. 1 mag at a distance of d fa 26 Mpc. SN 2004gq was 
also originally classified as Type Ic, but was later re-classified 
as Type lb after the appearance of strong helium features in 
a spectrum obtained on 2005 Jan 7.3 UT (Modjaz & Falco 
2005). As in the case of SN 2004dk, late-time spectroscopy by 
Maeda et al. (2008) indicated an ordinary hydrogen-poor SN 
lb spectrum. Detailed optical photometry for SN 2004gq was 
presented by Drout et al. (201 1) revealing a peak magnitude of 
M v fa -17.6 mag. 

We initiated VLA observations of SN 2004gq on 2004 Dec 
16 UT, roughly Af fa 8 days after explosion. Our observations 
of SN 2004gq are summarized in Table 4 spanning frequencies, 
v = 1.4,4.9,8.5 and 15 GHz. At 8.5 GHz, the SN reaches a 
peak luminosity of L v fa 5.5 x 10 27 erg s" 1 Hz" 1 at Af fa 21 
days. At Af fa 400 days, our observations reveal evidence for 
an abrupt fading at all frequencies. There are no subsequent 
radio observations available to determine if the emission later 
re-brightened. We estimate the abruptness factor of the radio 
steepening to be 0.21 and the power law decay to be at least as 
steep as F v oc f" 5,9 . 

3. A MODEL FOR THE RADIO EMISSION 

Supernova observations are typically limited to the optical 
and near-IR bands where the signal is dominated by thermal 
emission from the inner, slow-moving layers of ejecta that carry 
the bulk of the kinetic energy. On the other hand, radio ob- 
servations trace the outer, fastest-moving layer, the "shock- 



wave", to which only 0.01-0.1% of the total energy is coupled 
(e.g., Matzner & McKee 1999). Long duration GRBs are the 
notable exception, for which the relativistic jets and slower SN 
ejecta carry comparable energy, each approaching, and some- 
times exceeding, E ~ 10 51 erg (e.g., GRB 030329; Berger et al. 
2003b; Mazzali et al. 2003; Frail et al. 2005). 

As the Shockwave plows through the circumstellar medium, 
it shock-accelerates electrons into a power-law distribution of 
relativistic energies, N(-f) oc -f~ p , above a minimum Lorentz 
factor, 7 m . The Shockwave also generates amplified magnetic 
fields which cause the electrons to gyrate, emitting synchrotron 
radiation that peaks in the cm-band in the first days to weeks 
following the SN explosion. The synchrotron emission is sup- 
pressed at low frequencies due to absorption processes. In the 
case of SNe Ibc, Chevalier (1998) showed that the dominant ab- 
sorption process is internal synchrotron self-absorption (SSA), 
whereas external absorption processes (e.g. free-free) may con- 
tribute significantly for some SNe II plowing into denser cir- 
cumstellar environments with slower shock speeds (e.g., Frans- 
son & Bjornsson 1998). Self-absorption defines the spectral 
peak frequency for the synchrotron spectrum, v p , below which 
the radio spectral index is (3 = 5/2 (optically-thick) and above 
v p the spectral index is j3 = —(p- l)/2 (optically-thin; Rybicki 
& Lightman 1979). For SNe Ibc, a value of p fa 3 is typical 
(Chevalier 1998). As the Shockwave expands, the optical depth 
to SSA decreases, causing v p to cascade to lower frequencies 
with time. 

3.1. Radio Spectral Energy Distributions 

In our preliminary analysis of the SN data, we considered 
the multi-frequency radio observations from each epoch to ex- 
amine the properties of the synchrotron spectrum at various 
stages in the radio light-curve evolution. As shown by Cheva- 
lier (1998), a radio SN spectral energy distribution character- 
ized by SSA is well described as, F v = 1.582/v ( /f//f/ / ,) 5 / 2 [l - 

exp(-(^/f /J )" (/ ' +4)/2 )]. For each radio spectrum, we used a \ 2 
minimization technique to determine p, v p , and the associated 
flux density at the spectral peak, F VjP . Trie fitting routine was 
carried out independently at each epoch for which we obtained 
observations in at least three frequencies. The radio emission 
from all three SNe is well described by a SSA spectrum across 
multiple epochs (Figure 6). In each case, the electron energy 
index is consistent with p fa 3 and the optically thick spectral in- 
dex is not steeper than (3 fa 5/2 indicating that any contribution 
from external absorption processes (e.g., free-free absorption) 
is negligible. As can be seen from the Figures, v p cascades to 
lower frequencies with time and is accompanied by a decrease 
in F VtP . 

We make the reasonable assumption that the post-shock en- 
ergy density is equally fractioned between relativistic electrons 
(e e ) and amplified magnetic fields (e B ), i.e. equipartition, with 
the remaining energy coupled to shocked protons (e p ). In this 
framework, the radius of the emitting material may be esti- 
mated from the spectral parameters measured at each epoch, 
R«5.7x 10 15 (L„. ; ,/10 27 erg s" 1 YbT l fl l9 {v p /%.5 GHz)" 1 cm. 
Here we have assumed that the emitting region occupies half of 
the total volume enclosed within a spherical Shockwave. Com- 
bined with an estimate of the explosion date, the Shockwave 
radius enables an estimate of the time-averaged shock velocity, 
v. 
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Similarly, the strength of the amplified magnetic field may 
be estimated from the properties of the radio spectrum as B w 
0.93 (L^/10 27 erg s" 1 Hz" 1 )" 2 ' 1 V P /8.5 GHz) G. From this 
estimate, the total internal energy density of the radio emit- 
ting material is, U « B 2 /Sires- It is reasonable to assume that 
the partition fractions are each constrained to be (e e , cb) ^1/3 
(Chevalier 1998). A further constraint is implied by the shock 
velocity and the requirement that 7 m ;> 1 leads to a lower limit 

on the electron partition fraction, e e ;> 0.1(v/0.2c)~ 2 for p w 3 
(Soderberg et al. 2005; Chevalier & Fransson 2006; Soder- 
berg et al. 2010b). In line with the typical velocities measured 
for SNe Ibc Shockwaves, v « 0.2c, we make the reasonable as- 
sumption that e e = cb = 0.1 throughout this paper. 

3.2. A Dynamical Model 

We adopt the formalism of Soderberg et al. (2005) in model- 
ing the dynamics of the Shockwaves. In this model, the shock- 
wave parameters evolve self-similarly with time, R oc t a " and 
B oc t a ". As shown by Chevalier (1982), the expansion of 
the slightly decelerated Shockwave is characterized by a R = 
(n-3)/(n-s) where n describes the density profile of the ejecta, 
Pej oc r~ n , and s describes the density profile of the circumstel- 
lar medium, pcsm oc r~ s . The self-similar model is constrained 
such that n > 5 and in the case of a stellar wind density profile, 
5 = 2 and Ob = — 1- Using this general dynamical prescription, 
we fit the flux density measurements at each epoch and fre- 
quency simultaneously. 

The flux density is parametrized as 

/ \ (4a„-a B )/2 

in erg/cm 2 /s/Hz where F^ and Fj are Bessel functions and 
x = 2/3{v/v m ) with v m defining the critical synchrotron fre- 
quency, v m = 2>"f} n qB/(Ainn e c). The optical depth to SSA, t v , 
evolves with time as 

, t x Q+p/2)a B -«2p-3)a R -2(p-2) 

T ^ V)=C AT Q ) (2) 

x *T ( " +4)/2 F 2 (x) 

where C/ and C T are normalization constants of the flux den- 
sity and optical depth, respectively, and to is a reference epoch. 
By fitting for C/, C T , a R , and as, we derived estimates for the 
physical parameters, Rq, Bq, and j m fl at time, fo. 

The number density of emitting electrons may be estimated 
from these parameters as, n e « B 2 / {A-K^ m m e c 2 ) where we main- 
tain the assumption that p « 3 and e e = eg. Thus, the num- 
ber density of shocked electrons scales temporally as a„ e = 
2aB-ay and radially as s = -(2cxb~ a 7 )/a R . The density of 
unshocked electrons in the circumstellar medium is related by 
the shock compression factor, rj, such that «csm = n e /rj. The 
progenitor mass loss rate follows directly as M ss 2im e m p R 2 v w 
where we have assumed a compression factor of r\ = 4 and a 
helium-rich stellar wind (nucleon-to-electron ratio of 2). Here, 
v w is the velocity of the progenitor wind and we adopt v w = 
10 3 km s" 1 as measured for Galactic Wolf-Rayet stars (Cappa, 
Goss & van der Hucht 2004; Crowther 2007). The temporal 



evolution of the inferred mass loss rate is, M oc t 2a " +2 , and is 
constant for a standard stellar wind profile with a steady wind 
velocity (s = 2). 

3.3. Modeling Results 

For each SN, we model the early radio data only, excluding 
the unusual late-time light-curve behavior in the fit. The tem- 
poral and spectral evolution of the SNe are well described by 
the dynamical model (Figures 2-5). For each SN, we report the 
two model parametrizations which best represent the data (see 
Table 5). In the primary fit, we used a \ 2 minimization tech- 
nique to determine reasonable values for C/, C T , a R , and Og. In 
the "standard" model fit, we fixed s = 2, olb = -1 and a R = 0.9 
as expected in the typical scenario where the forward shock is 
only slightly decelerated (Chevalier & Fransson 2006), and fit 
only for C/ and C T . We extracted estimates for the physical 
parameters (R, B, E, and n e ) and their associated temporal and 
radial evolution. The results are detailed below for each SN 
and summarized in Table 5. We compare the temporal and ra- 
dial evolution of the physical parameters for these three SNe in 
Figure 7. 

In the case of SN2004cc, the best fit is characterized by 
cyr « 0.9, and a steep evolution of the magnetic field is required, 
as ~ — 1.4. Together these, imply a steep CSM density gradi- 
ent, s ~ 2.8. At to = 10 days, both models predict a velocity 
of vo ~ 0.1c. However, the models differ regarding the energy, 
with the best fit model predicting an energy of £b ~ 1 -0 x 10 48 
erg while for the standard model, £b ~ 2.7 x 10 47 erg. Assum- 
ing the wind velocity is unchanged, the inferred mass loss rate 
varies with time; we estimate Mq ss 1.3 x 10~ 4 M© yr" 1 at to and 
just prior to the late-time re-brightening episode it is a factor of 
- 10 lower, M = 2.1 x 10" 5 M© yr -1 . 

For SN 2004dk we find that the standard model provides an 
excellent description of the early radio data; our primary fit is 
not dissimilar from the standard model. We estimate a time- 
averaged Shockwave velocity of vo ~ 0.2c and an energy of 
£0 ~ 1-7 x 10 47 erg at to = 10 days. We find the progenitor 
mass loss rate to be constant, Mo ~6.0x 10~ 6 M© yr" 1 . 

In the case of 2004gq, we find that the best fit is nearly con- 
sistent with a standard model, but with a shallower density pro- 
file, s = 1.5. In both models, the physical parameters at fo = 10 
days are essentially the same. The implied time-averaged ve- 
locity is, vo w 0.25c and the energy is £0 ~ 5 x 10 47 erg. The 
implied progenitor mass loss rate is M ~9x 10~ 6 M© yr -1 
and is not constant in the best fit model. Just prior to the 
epoch of the late-time flux steepening, the mass loss rate is 
M«4x 10" 5 M© yr" 1 . 

4. LATE-TIME FLUX VARIATIONS 

Intriguingly, the late-time data for all three SNe reveal un- 
usual behavior, in the form of an strong radio re-brightening for 
SNe 2004cc and 2004dk, and an abrupt fading for SN 2004gq. 
These late-time data were not included in the model fitting rou- 
tines described above. For each SN, we measure the flux varia- 
tion by comparing the late-time measurements (or limits) to an 
extrapolation of the model. In Table 5 we report the late-time 
flux modulation factor for each SN. 

Radio variability in young SNe is typically caused by mod- 
ulations in the circumstellar density profile, but engine-driven 
SNe may also produce radio variability due to off-axis ejecta 
components (jets) or prolonged central engine activity (see 
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Soderberg et al. 2006a for a discussion of radio SN variability). 
In the case of CSM density modulations, the flux variations are 
accompanied by a significant and abrupt shift in the location 
of synchrotron self-absorption frequency. This unique obser- 
vational clue has enabled circumstellar density variations to be 
inferred for several core-collapse SNe exhibiting radio variabil- 
ity (Ryder et al. 2004; Soderberg et al. 2006a). 

Chevalier (1998) show that for SNe dominated by SSA, 
F v oc R 2 NqB (p+v> / 2 , where JVo is the normalization of the elec- 
tron energy distribution at j m « 1. Assuming eg is unchanged 
following an abrupt modulation in the circumstellar density, we 
have B 2 oc n e v 2 and thus expect roughly F v oc ni 1 * 5 ^ 4 or F„ oc n 2 
for p i=a 3. Doubling the CSM density therefore produces a fac- 
tor of ~ 4 increase in the optically-thin flux density. Assuming 
that the wind velocity is unchanged (a reasonable assumption 
since v w is roughly set by the escape velocity) and the CSM 
density profile is characterized by s = 2, the implied mass loss 
rate modulation scales similarly, AF V oc AM 2 . In cases where 
s 7^2, we adopt AF V oc A(M) 2 where (M) represents a radially- 
averaged mass loss rate. 

4. 1 . Mass Loss Variations 

SN2004cc shows evidence for a late-time re-brightening, 
nearly 100 times the flux density expected from an extrapola- 
tion of the model fits to the early data (Figure 2). The spectral 
index between frequencies of v = 8.5 and 15 GHz at At « 130 
days is j3 « -0.1, indicating that the spectral peak was be- 
tween these two frequencies at that time. However, v p cas- 
caded through both frequencies several months earlier, indicat- 
ing that the re-brightening is associated with an increase of v p . 
We therefore attribute the radio re-brightening to a significant 
circumstellar density modulation. A factor of <~ 1 1 increase 
in the optically-thin flux corresponds to a CSM density jump 
by a factor of ~ 3.3 and a mass loss rate of (M) « 6.9 x 10" 5 
M Q yr" 1 . 

More than four years after explosion, an abrupt re- 
brightening (factor of ~ 40) is observed for SN 2004dk and 
may be associated with a flattening of the spectral index be- 
tween v = 4.9 and 8.5 GHz, 8 « -0.2 (Fig. 4) 6 . Similar to the 
case of SN 2004cc, the light-curves at these frequencies peaked 
much earlier and hence the re-brightening is associated with an 
upward shift of v p . We attribute this re-brightening to an in- 
crease in the circumstellar density by a factor of <~ 6. A density 
modulation of this size implies an increased mass loss rate of 
(M}^3.5xl0- 5 M e yr 1 . 

In the case of SN 2004gq, the late-time observations indicate 
a dramatic decline in the light-curves at Af « 400 days. At this 
point, the flux density was below our detection threshold so a 
spectral index measurement was not possible, however, we note 
that the steepness of the decline (F v oc f" 6 ) is roughly consistent 
with the evolution expected for spherical adiabatic expansion 
losses of particles with an energy distribution of p 3 and flux 
freezing (Shklovskii 1960; Kellermann & Pauliny-Toth 1968). 
Attributing the flux decline to a density drop, the measured up- 
per limits imply a drop in the CSM density by a factor of ;> 5. 
We constrain mass loss rate associated with this lower density 
wind to (M) £ 1 .8 x 10" 5 M Q yr" 1 . 

In Figure 7, we compare the circumstellar density profiles for 

6 We note that the S/N of the late-time observations makes it difficult to precisely 
sensitivity of the EVLA are planned to follow the long-term evolution of the radio 



all three SNe as inferred from our dynamical modeling and in- 
cluding the abrupt modulations. Our radio observations point to 
significant (factor of ~ 3 - 6) circumstellar density modulations 
on radial scales of r « (1-50) x 10 16 cm. Assuming a typical 
progenitor wind velocity of v w = 10 3 km s" 1 , these density mod- 
ulations are traced back to mass loss variations spanning a few 
years to a century prior to explosion. 

4.2. A Comparison to Other Radio SNe 

Our modeling of the early radio light-curves for SNe 2004cc, 
2004dk, and 2004gq indicate mass loss rates spanning M « 
(0.6- 13) x 10" 5 M yr" 1 . However, these estimates vary by 
factors of ~ 3 - 6 following the observed late-time flux den- 
sity modulations. In comparison with other radio SNe simi- 
larly observed on sub-parsec scales, these SNe represent exam- 
ples of the strongest flux modulations ever reported. To date, 
such large-scale (factor of ^5-10) radio modulations have 
only been reported for a few SNe, including SN 1987A (Ball 
1995; Zanardo et al. 2010), SN 1996cr (Baueref al. 2008), 
SN 2001 em (Stockdale et al. 2005), and SN2007bg (Soder- 
berg 2009). 

Small-scale mass loss variations (factor of <~ 2) are more 
common and have been inferred for several core-collapse SNe 
including 1979C (Type IIL; Weiler et al. 1991), 200 lig (Type 
lib; Ryder et al. 2004), 2003bg (broad-lined Type lib; Soder- 
berg et al. 2006a), and 2008ax (compact progenitor Type 
Ib/IIb; Roming et al. 2009; Chevalier & Soderberg 2010). Ra- 
dio variability was also clearly seen for the relativistic central 
engine-powered 1998bw (associated with GRB 980425; Kulka- 
rni et al. 1998) and 2009bb (Soderberg et al. 2010a; Bi- 
etenholz et al. 2010), both of which are broad-lined SNe Ic. 
However, it remains debated whether the variability was pro- 
duced by central engine activity, additional ejecta components, 
or variations in the CSM density profile. Therefore, radio vari- 
ability extends across the various SN Ibc sub-classes, includ- 
ing helium-rich, helium-poor, broad-lined, engine-driven and 
even hydrogen-rich (Type cllb; Chevalier & Soderberg 2010). 
We stress, however, that density-induced flux variations are 
not ubiquitous as evidenced by several well-studied SNe with 
smooth radio light-curves including SN 1993J (lib; (Weiler et 
al. 2007)), SN2003L (Ic; Soderberg et al. 2005), SN2008D 
(lb; Soderberg et al. 2008), and GRB-SN 2003dh (i.e. the radio 
afterglow of GRB 030329; Berger et al. 2003b), each observed 
on similar sub-parsec radial scales. 

We compare the radio-derived mass loss rates for SNe 
2004cc, 2004dk, and 2004gq with those of other local SNe Ibc 
and engine-driven (GRB-SN) explosions within a comparable 
volume, d <; 150 Mpc. The comparison sample of SNe repre- 
sents a literature compilation (see Figure 1 caption), with in- 
ferred mass loss rates spanning M « 10~ 7 - 10~ 2 M Q yr" 1 and 
a distribution peak nearM « 10~ 5 M Q yr" 1 (see also Chevalier 
& Fransson 2006). We derive mass loss rates for SNe 2004cc, 
2004dk, and 2004gq that overlap with the higher end of this 
SN Ibc distribution. In contrast, the engine-driven explosions 
considered here populate the low end of the mass loss rate dis- 
tribution, M ss (1-50) x 10" 7 M yr" 1 . 

A modest dispersion in the M distribution may be expected 
due to variations in the shock partition fractions (e e , eb), how- 

/ measure the spectral index. Further observations of SN 20()4dk with the increased 
emission. 
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ever, these are unlikely to deviate strongly (e.g., J> 100) from 
equipartition (see Chevalier & Fransson (2006) for discussion). 
Indeed, in the case of SN2011dh the partition fractions were 
found to have a ratio e e /e B ~ 30 (Soderberg et al. 2012). For 
deviations from equipartition, the mass loss rates reported here 
would scale as M oc e^^^eg 11 ^ 19 . We emphasize, however, that 
this only increases the inferred mass loss rates. 

Some dispersion in the M distribution is also expected from 
an intrinsic spread in progenitor wind speeds; we emphasize 
that the radio observations only constrain the quantity, (M/v w ), 
and we have assumed a constant wind, v w = 10 3 km s" 1 for each 
SN. However, v w is regulated by the progenitor escape speed 
and is not seen to vary by more than a factor of a few in the set 
of well-studied Galactic WRs (see §5). Therefore, the factor 
of 10 5 spread in the inferred late mass loss rates (last 10-100 
years) of SN Ibc progenitors points to an intrinsically broad 
dispersion 7 . 

In comparison to other supernova sub-types, the mass loss 
rates for these SN Ibc progenitors are lower than those inferred 
for most Type Iln SNe, revealing intense progenitor mass loss 
episodes with M « 0.01 - 1 M Q yr" 1 and wind velocities of 
v w w 300 km s" 1 in the decades leading up to explosion (Smith 
2008; Kiewe et al. 2012). The local CSM densities of SNe Iln 
are so high that the radio emission typically becomes optically- 
thin only years after the explosion in the cm-band such that the 
early radio signal is suppressed below detectable limits (e.g., 
SN1988Z; Williams et al. 2002). 

5. THE NATURE OF CSM DENSITY MODULATIONS 

While it is generally understood that the massive progenitor 
stars of SNe Ibc (including relativistic SNe) lose their envelopes 
prior to explosion, the mechanism(s) by which the material is 
removed and the relevant timescale(s) remains debated (Smith 
et al. 201 1). Here we review the possible causes of circumstel- 
lar density modulations in core-collapse SNe. 

Isolated Wolf-Rayet stars, the favored progenitors of SNe 
Ibc, are seen to lose mass at a typical rate of M m (0.5-2) x 
10" 5 M yr" 1 (including Galactic WN, WC, and WO subtypes; 
Cappa, Goss & van der Hucht 2004; Crowther 2007) through 
line-driven winds traveling at v w w (0.7-2.5) x 10 3 km s" 1 , pro- 
pelled by radiation pressure and momentum coupling (Castor, 
Abbott & Klein 1975; Conti 1978). These winds shape the im- 
mediate circumstellar environment on sub-parsec scales. The 
strength of such line-driven winds scales almost linearly with 
metallicity (M oc Z° 8 ; Vink & de Koter 2005), and an upper 
limit on this process is set by saturation of the metal lines at 
M « 10" 4 M Q yr" 1 . Small-scale clumping within the winds due 
to turbulence can produce moderate density modulations (factor 
of 2-4; see Moffat 2008). At similar circumstellar radii, mod- 
est density modulations may be present due to subtle variations 
in the stellar wind velocity. The collision of winds in a tight 
binary systems can also give rise to CSM density modulations 
in the immediate vicinity of the progenitor. 

At larger circumstellar radii (a few parsec) an abrupt density 
jump is likely produced by a wind termination shock between 
the fast Wolf-Rayet wind and the slower red supergiant (RSG) 
wind (Garcia-Segura, Langer & Mac Low 1996). Numerical 

7 We note that the subset of rare, radio luminous SNe are likely over-represented 
Soderberg et al. 2005; Chevalier & Fransson 2006). A distance-limited study of rac 
the intrinsic dispersion in M. This is the focus of a separate paper. 



simulations indicate that the structure of the wind-termination 
shock is shell-like, and that the interaction of the SN shock- 
wave with this density feature gives rise to a strong radio and 
X-ray signal, accompanied by hydrogen recombination lines in 
the optical band (e.g., Dwarkadas, Dewey & Bauer 2010). 

Massive circumstellar shells are also produced in binary pro- 
genitor systems where a common envelope is ejected by the 
motion of a companion star (Podsiadlowski, Joss & Hsu 1992). 
The envelope material moves more slowly (v w « 10 km s" 1 ) 
than the Wolf-Rayet wind and so becomes accelerated and com- 
pressed into a thin dense shell. The interaction of the shock- 
wave with the common envelope shell produces an abrupt in- 
crease in the synchrotron emission and Balmer series recombi- 
nation lines visible in optical spectra. Late-time radio and op- 
tical observations of SN 2001 em revealed these characteristics 
(Soderberg, Gal- Yam & Kulkarni 2004; Stockdale et al. 2005), 
motivating Chugai & Chevalier (2006) to suggest that a com- 
mon envelope was lost at a rate of (2-10) x 10" 3 M Q yr" 1 and 
located at r « 7 x 10 16 cm at the time of the explosion. 

Large density modulations may also be fueled by violent stel- 
lar outbursts such as those observed from Luminous Blue Vari- 
able stars (LBVs; Humphreys & Davidson 1994). LBVs are 
thought to represent a short-lived phase in the evolution of some 
supergiants, and immediately precede the development of clas- 
sic WR features (Crowther 2007). Along this line, LBVs and 
WRs both show evidence for strong stellar winds and H-poor 
envelopes. LBVs give rise to giant mass loss ejections exceed- 
ing the saturation limit associated with line-driven winds and 
are therefore attributed to metallicity-independent continuum- 
driven winds and/or hydrodynamic eruptions (Smith & Owocki 
2006). The outbursts result in dense and massive circumstel- 
lar shells that expand into a strong (and quasi-steady) wind. 
Multi-wavelength emission is observed from the resulting neb- 
ulae surrounding local LBVs on radial scales of of r <; 1 pc 
(e.g., Umana et al. 2010). 

The long-lived and luminous multi-wavelength emission 
from Type Iln SNe is commonly attributed to the abrupt in- 
teraction of the Shockwave with LBV-like mass ejections in the 
decades leading up to the explosion (e.g., Chugai et al. 2004). 
A recent (and extreme) example is the Type Iln SN 2006gy, for 
which a progenitor mass loss rate of M w 0.5 M Q has been in- 
ferred (Ofek et al. 2007; Smith et al. 2007; see van Marie 
et al. 2010 for a review). An intriguing aspect of the SN Iln 
mass loss histories is the synchronization required between the 
violent ejection episodes and the ultimate explosion (see Smith 
2008 for a discussion). An LBV-like eruption seen from the 
Type lb SN 2006jc progenitor star two years prior to explosion 
similarly indicates a synchronization, and in turn, a direct link 
between SN Ibc and Iln progenitors (Foley et al. 2007; Pas- 
torello et al. 2007). At longer wavelengths, a link was also 
suggested by attributing the episodic radio variability observed 
for Type Ib/IIb SNe 2001ig and 2003bg to S-Doradus pulsations 
of the progenitor star (Kotak & Vink 2006, but see Soderberg 
et al. 2006a and Chevalier & Soderberg 2010 for a different 
interpretation). 

5.1. SNe 2004cc, 2004dk, 2004 gq 

in this literature compilation since they are preferentially detected (e.g., SN 2003L; 
Jio SN luminosities, and in turn, progenitor mass loss rates is required to determine 
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We consider the observed multi-wavelength properties for 
SNe 2004cc, 2004dk, and 2004gq in light of the various causes 
of circumstellar density modulations discussed above. As dis- 
cussed in §4, the late-time radio data for these three SNe imply 
mass loss rates spanning M ss (2-7) x 10~ 5 M yr" 1 at cir- 
cumstellar radii of r « (1 -50) x 10 16 cm. Their late-time radio 
variability implies circumstellar density modulations by a fac- 
tor of ~ 3-6. Moreover, early optical photometry and late-time 
(nebular) spectroscopy are consistent with observations of typi- 
cal SNe Ibc and do not reveal any evidence for unusually strong 
circumstellar interaction (e.g., H recombination lines; §2). 

The intermediate-level density modulation factors and ordi- 
nary optical properties observed for these SNe indicate that it 
is unlikely that their Shockwaves are interacting with a H-rich 
common envelope. Attributing the abrupt density modulations 
to a clumpy wind would imply that the clumps are exceed- 
ingly large (comparable to the size of the expanding shock- 
wave) in order to produce an observable upward shift in the 
self-absorption frequency. This is inconsistent with the ba- 
sic picture for small-scale, turbulence-driven WR wind clumps 
(e.g. Moffat 2008). Interaction with the wind termination 
shock should occur at significantly larger radii (several pc) so 
we find this scenario unlikely. However, it is possible to reduce 
the radius of the termination shock if these SNe are embedded 
in highly pressurized environments (Chevalier, Li & Fransson 
2004) or if their progenitors were runaway stars prior to explo- 
sion (Cantiello et al. 2007; van Marie et al. 2006). 

In comparison to mass loss rates measured for Galactic Wolf- 
Rayet stars, the inferred circumstellar densities inferred for 
these three SNe are at the high end of the distribution (Fig- 
ure 8). They show overlap with the mass loss rates inferred for 
radio luminous SNe 2003L and 2003bg on similar radial scales 
(M sa 1 - 3 x 10" 4 M yr -1 ; Chevalier & Fransson 2006), which 
are roughly consistent with the saturation limit for line-driven 
winds. 

This, together with the abruptness of the mass loss varia- 
tions and their apparent synchronization with respect to the 
explosion date, resembles the properties of SNe Iln. Their 
unusual mass loss properties may indicate some contribution 
from metallicity-independent envelope-stripping mechanisms 
(e.g., continuum-driven winds and/or hydrodynamic eruptions). 
While the density enhancements inferred here are not as ex- 
treme as those seen for Galactic LBVs or inferred for SN Iln 
progenitors, their abrupt nature is reminiscent. We note that 
recent theoretical efforts suggest that such abrupt mass loss 
episodes could be driven by pulsations or gravity waves induced 
by convective motions (Yoon & Cantiello 2010; Quataert & Sh- 
iode 2012). 

5.2. A Metallicity Dependence? 

A key question raised by these observations is whether 
the circumstellar environments inferred for some SNe are 
shaped by metallicity-dependent line-driven winds or "LBV- 
like" metallicity-independent mass loss mechanisms. While it 
is not possible to measure the metallicity of the progenitor star 
post-explosion, the metallicity of the explosion site on a sub- 
kpc scale may serve as a reasonable proxy for that of the short- 
lived progenitor system. A comparison between the explosion 
site metallicities and the radio-derived mass loss rates tests the 



metallicity-dependence of the envelope-stripping mechanism. 

As a rough estimate of the explosion site metallicity, we 
adopt the absolute B-band magnitudes for the three host 
galaxies (Table 1) and the luminosity-metallicity relation of 
Tremonti et al. (2004). We find host galaxy metallicities of 
[12+log(0/H)]py 8.9, 8.8, and 9.0 for SNe 2004cc, 2004dk, and 
2004gq, respectively. These estimates are somewhat higher 
than the solar metallicity, [12+log(0/H)]w 8.69 (Allende Pri- 
eto, Lambert & Asplund 2001). However, these values are 
global metallicity estimates, and there is evidence that the ex- 
plosion site metallicities of SNe Ibc can be significantly differ- 
ent than those estimated globally due to metallicity gradients 
within the host galaxy (see Anderson et al. 2010; Modjaz et al. 
201 1 for a discussion). Along this line, Anderson et al. (2010) 
recently reported an explosion site metallicity for SN 2004gq 
of [12+log(0/H)]w 8.8, lower than that estimated globally for 
the host galaxy 8 . Explosion site metallicities are currently un- 
available for SN2004cc and SN2004dk; however, Modjaz et 
al. (201 1) estimate that the offset between global and sub-kpc 
metallicities for SNe Ibc is typically of order ±0.2 dex. We con- 
clude that the explosion sites metallicities for these three SNe 
are likely typical, probably not far from the Solar value and 
consistent with the range of metallicities at which line-driven 
winds operate efficiently. 

Finally we note that the current range of explosion site metal- 
licities spans just 0.5 dex for ordinary SNe Ibc and engine- 
driven explosions (Modjaz et al. 2011). However, the dis- 
persion in radio-derived mass loss rates spans several orders 
of magnitude, significantly broader than expected if the M-Z 
relation is roughly linear. This preliminary comparison sug- 
gests that metallicity-dependent mass loss processes (e.g., line 
driven-winds) are not the only mechanism by which the pro- 
genitor envelope is removed in SNe Ibc and/or engine-driven 
explosions. 

6. CONCLUSIONS 

We present extensive radio observations and detailed mod- 
eling for the three nearby Type Ibc SNe 2004cc, 2004dk, and 
2004gq. We show that the physical properties of their shock- 
waves are typical of other SNe Ibc while their circumstellar en- 
vironments are characterized by unusually strong density mod- 
ulations. We conclude with the following points: 

• The radio spectra for all three SNe are well described by 
a self-absorbed synchrotron spectrum with a relativistic 
e" population characterized by 7,,, ~ 1 -2 and p « 3. 

• Detailed modeling of the early radio data shows them to 
be consistent with free-expansion model in which a non- 
relativistic shock interacts with a stellar wind environ- 
ment. We extract Shockwave velocities of v <~ 0. 1 -0.25c 
and energies of E w (2-10) x 10 47 erg at At « 10 days. 

• Each SN is characterized by a late-time radio re- 
brightening or fading which we attribute to abrupt and 
significant (factor of <~ 3 - 6) density modulations at cir- 
cumstellar radii of, r « (1 -50) x 10 16 cm. The in- 
ferred mass loss rates rates lie at the high end of the M- 
distribution measured for Galactic WRs 



8 We note, however, that the explosion site metallicity was estimated using the emission line metallicity diagnostic of Pettini & Pagel (2004) which prevents a 
direct comparison of the two Z-estimates. 
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• If we attribute the density modulations to variations in 
the progenitor mass loss rates then they approach the sat- 
uration limit for line-driven winds. If the density modu- 
lations are instead due to the collision with a wind termi- 
nation shock, the SNe must reside in highly pressurized 
environments or be runaways in order to squeeze the ter- 
mination shock to small circumstellar radii. 

• The possible synchronization of the mass loss variations 
with the SN explosion is reminiscent of SN Iln progeni- 
tors. 

• A comparison of the radio-derived M estimates and ex- 
plosion site metallicities for stripped envelope SNe could 
shed light on the primary mass loss mechanism(s). 

Finally, we note that that the strong flux density variations 
reported here are unlikely to be ubiquitous since they are recov- 
erable by radio campaigns designed to search for off-axis GRB 
jets associated with SNe Ibc on timescales of months to years 



after explosion (Perna & Loeb 1998; Paczynski 2001). Such ra- 
dio searches have uncovered strong late-time variability in less 
than ~ 10% of all SNe Ibc (Berger et al. 2003a; Soderberg et 
al. 2006c). Future radio efforts will shed light on the variabil- 
ity of radio SNe over various time baselines and therefore able 
to probe density modulations at a wide-range of circumstellar 
radii. In particular, the improved sensitivity of the Expanded 
Very Large Array (Perley et al. 2009) will enable detailed mon- 
itoring of radio SNe to longer timescales, enabling the search 
for variable progenitor mass loss histories at larger circumstel- 
lar radii, possibly associated with earlier pre-explosion epochs. 
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for Undergraduates (REU) and Department of Defense Awards 
to Stimulate and Support Undergraduate Research Experiences 
(ASSURE) programs under Grant no. 0754568 and by the 
Smithsonian Institution. 
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Fig. l. — Radio light curves for 2004cc (blue), 2004dk (green), 2004gq (red), are compared with other Type Ibc supernovae, 
including SNe 1983N (Sramek, Panagia & Weiler 1984), 1990B (van Dyk et al. 1993), 19941 (Weiler 2006), 2001ig (Ryder et 
al. 2004), 2002ap (Berger, Kulkarni & Chevalier 2002), 2003L (Soderberg et al. 2005), 2003bg (Soderberg et al. 2006a), 2007gr 
(Soderberg et al. 2010b), 2008D (Soderberg et al. 2008), 2008ax (Roming et al. 2009), and 201 ldh (Krauss et al. 2012) all 
shown in grey. The radio afterglows of nearby engine-driven explosions within a similar volume are shown in black, including 
GRB-SN 1998bw (circles; Kulkarni et al. 1998), XRF-SN2006aj (squares; Soderberg et al. 2006b), and SN2009bb (triangles; 
Soderberg et al. 2010a). SNe 2004cc, 2004dk, and 2004gq are similar in radio luminosity to ordinary radio SNe and less luminous 
than engine-driven explosions on a similar timescale. 
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Fig. 2. — Radio observations of SN2004cc are shown various frequencies and compared to our dynamical models. Dashed lines 
represent the model parametrization assuming a density profile of n e oc r" 2 . We find a significantly better fit (solid lines) associated 
with a steeper density profile. Encircled data points indicate an unusual late-time re-bightening and were not included in the model 
fits. 
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Fig. 3.— Radio observations of SN 2004dk are shown various frequencies and compared to our dynamical models. The solid lines 
represent the best fit model parametrization, while the dashed lines represent a model in which the density profile to fixed as n e oc r~ 2 . 
Encircled data points at indicate an unusual late-time re-bightening and were not included in the model fits, shown in detail in Fig. 4. 
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Fig. 4.— Radio observations of SN 2004dk at late time, during the rebrightening episode, are shown at several frequencies. The SN 
is optically thin throughout the rebrightening. 
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Fig. 5.— Radio observations of SN 2004dk are shown various frequencies and compared to our dynamical models. The solid lines 
represent the best fit model parametrization, while the dashed lines represent a model parametrization in which the density profile is 
fixed to n e oc r~ 2 . The deep limits at late-time indicate an unusual steepening of the radio flux and were not included in the model fits. 
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Fig. 6. — Multi-frequency radio observations of SNe 2004cc, 2004dk, and 2004gq are well described by an SSA spectrum with 
F v oc v 5 l 2 in the optically thick regime and F v oc v~ l in the optically thin regime, associated with an electron index of p m 3. 
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Fig. 7. — The derived temporal and radial properties of the Shockwave and local environment are compared for SNe 2004cc (blue), 
2004dk (green), and 2004gq (red). The solid lines represent the best-fit parametrization of our dynamical model while the dashed 
lines correspond to a parameterization with a CSM denisty profile fixed to n oc r" 2 . 
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Fig. 8. — The density profile of shocked particles is shown as a function of circumstellar radius for SNe 2004cc (blue), 2004dk 
(green), and 2004gq (red). Modeling of the early data imply density values shown by the color dotted lines. The abrupt deviations 
from these fits at later epochs are associated with CSM density modulations. Radially averaged mass loss rates are shown for 
comparison (black dashed lines). The grey hatched region marks the observed range of mass loss rates for Galactic Wolf-Rayet stars. 
Assuming a wind velocity of 1000 km/s, these density modulations could be attributed to variable progenitor mass loss on a timescale 
of ~ 10- 100 years before explosion (top axis). 
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Table 1 
SN Sample Properties 



Name 


Host Galaxy 


Distance 1 


Host Magnitude 


Disc. Date 2 


Explosion Date 


Type 


Fpfi.5 GHz 


tp,&.5 GHz 






(Mpc) 


(M B ) 


(UT) 


(UT) 




(mJy) 


(days) 


2004cc 


NGC 4568 


18±2 


-19.6±0.3 


2004 Jun 10 


2004 Jun 2 


Ic 


4.5 


32 


2004dk 


NGC6118 


23 ±2 


-19.4 ±0.3 


2004 Aug 1 


2004 Jul 30 


lb 


2.5 


14 


2004gq 


NGC 1832 


26±5 


-20.1 ±0.7 


2004 Dec 1 1 


2004 Dec 8 


lb 


6.9 


21 



'We adopt cosmology independent distances from NED- ID when available, otherwise we adopt the model-derived NED distances 
assuming Hq = 73 km s" 1 Mpc" 1 , Q,m = 0.27, and 17a = 0.73. 

2 We estimate the explosion date by the linear mean of the UT dates bridging the optical SN discovery and the most recent pre- 
discovery non-detection. 



Table 2 
Observations of SN 2004cc 



Date 


Fv,%.5 GHz 


F u ,\5 GHz 


Fv.12.5 GHz 


F V ^MA GHz 


VLA 


(UT) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


Config. 


2004 Jun 17.2 






7.35 ±0.18 




D 


2004 Jun 18.1 


0.47 ±0.12 




8.33 ±0.16 




D 


2004 Jun 19.0 


0.93 ± 0.07 




9.44 ±0.11 


11.79 ±0.55 


D 


2004 Jun 24.2 


2.44 ±0.10 




13.83 ±0.29 


9.60 ± 0.57 


D 


2004 Jun 26.1 


3.66 ±0.10 


10.69 ± 0.37 






D 


2004 Jun 28.1 


3.87 ± 0.09 




11.29 ±0.30 


7.70 ±0.61 


D 


2004 Jun 30.1 


4.60 ± 0.09 


8.44 ± 0.21 






D 


2004 Jul 4.0 


4.13 ±0.08 




4.81 ±0.24 


1.89 ±0.63 


D 


2004 Jul 16.0 


<2.36 


1.36 ±0.40 


< 0.31 




D 


2004 Jul 20.0 


< 1.56 


0.64 ±0.15 


<0.38 


<0.61 


D 


2004 Aug 1.9 


0.48 ± 0.05 




0.72 ± 0.29 




D 


2004 Oct 14.6 


2.71 ±0.05 








A 


2004 Nov 21.6 


1.46 ±0.06 


1.37 ±0.18 


<0.21 




A 


2009 Feb 3.3 


<0.05 








B 
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Table 3 
Observations of SN 2004dk 



Date 


FvA.9 GHz 


Fv$.S GHz 


F u ,\5 GHz 


F v fi2.5 GHz 


VLA 


(UT) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


Coring. 


2004 Aug 7.1 




1.30 ±0.05 






D 


2004 Aug 11.0 


0.87 ± 0.06 


2.11 ±0.06 




2.78 ± 0.33 


D 


2004 Aug 13.0 


1.09 ±0.06 


2.42 ± 0.06 


3.57 ± 0.26 


2.72 ± 0.37 


D 


2004 Aug 18.lt 


0.95 ± 0.07 


1.76 ±0.07 


0.83 ± 0.23 


1.75 ±0.46 


D 


2004 Sep 2.0 


1.80 ±0.09 


1.56 ±0.13 


0.79 ±0.21 


0.69 ±0.18 


D 


2004 Sep 10.0 


1.06 ±0.07 


0.81 ±0.15 




<0.43 


A 


2004 Sep 18.0 


0.86 ±0.10 


1.19 ±0.11 




<0.47 


A 


2004 Oct 4.1 


0.46 ±0.15 


0.41 ± 0.09 


<0.34 


< 1.41 


A 


2004 Oct 17.8 


0.56 ±0.10 


0.47 ± 0.12 


<0.92 




A 


2005 Feb 12.4 


<0.26 


<0.17 






BnA 


2009 Feb 12.6* 


0.22 ± 0.06 








B 


2009 Feb 24.5* 


0.21 ±0.03 


0.19 ±0.03 






B 


2009 Apr 2.6 


0.34 ± 0.07 


<0.08 






B 


2009 Apr 19.0* 


0.21 ± 0.02 


0.15 ±0.02 






C 


2009 Sep 19.8 


0.28 ± 0.05 


0.17 ±0.02 




<0.57 


B 


2009 Oct 26.8 








<0.27 


D 



tPoor weather on 2004 Aug 18.1. 

*Flux density values reported by Stockdale et al. (2009). 

■^Observations also conducted at 1.4 GHz revealing a flux density of 0.18 ± 0.09 mJy. 
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Table 4 
Observations of SN 2004gq 



Date 


F V ,1A GHz 


^i/,4.9 GHz 


^V,8.5 GHz 


F u .15 GHz 


VLA 


(UT) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


Config. 


2004 Dec 16.2 






2.00 ± 0.05 




A 


2004 Dec 19.3 






3.00 ± 0.06 




A 


2004 Dec 24.3 




1.81 ±0.13 


5.61 ±0.10 




A 


2004 Dec 27.2 




2.06 ±0.14 


5.68 ±0.11 


6.27 ± 0.22 


A 


2005 Jan 8.2 




5.75 ± 0.09 


6.40 ± 0.06 


3.53 ± 0.20 


A 


2005 Jan 15.2 


<0.92 


4.90 ±0.10 


4.47 ± 0.05 


1.96 ±0.17 


A 


2005 Jan 21.1 


0.85 ± 0.24 


4.12 ±0.09 


3.43 ±0.10 




BnA 


2005 Jan 28.1 




4.45 ±0.10 


2.95 ± 0.08 


1.54 ±0.50 


BnA 


2005 Feb 8.2 


1.49 ±0.27 


3.56 ± 0.08 


2.40 ± 0.06 


1.64 ±0.15 


BnA 


2005 Feb 27.1 


2.62 ±0.13 


3.44 ±0.12 


1.91 ±0.07 




B 


2005 Mar 13.0 


2.45 ± 0.55 








B 


2005 Mar 20.0 


2.75 ±0.16 


2.68 ±0.12 


1.33 ±0.09 


0.27 ± 0.08 


B 


2005 Apr 13.0 




1.84 ±0.09 


1.15 ±0.07 


0.63 ±0.18 


B 


2005 May 21.9 


3.21 ± 0.29 


0.97 ± 0.09 


0.64 ± 0.06 


<0.66 


B 


2005 Jun 5.8 


2.18 ±0.17 


0.91 ±0.07 


0.59 ± 0.06 


<0.67 


B 


2005 Aug 7.6 


2.62 ± 0.30 


0.57 ± 0.30 


0.45 ± 0.07 


<0.55 


C 


2005 Aug 27.6 


2.09 ±0.31 


1.13 ±0.29 


0.41 ± 0.05 


<0.47 


C 


2005 Oct 1.5 


2.71 ± 0.20 


<0.62 


0.35 ± 0.04 




C 


2005 Oct 27.4 




<0.48 


0.43 ±0.14 




DnC 


2006 Jan 2.3 






0.41 ± 0.05 




D 


2006 Mar 26.9 


<0.10 


<0.09 


<0.14 


< 1.17 


A 


2009 Apr 6.0 






<0.08 




B 



23 



Table 5 
Model Parameter Fits 





2004cc 




2004dk 




2004gq 






Best Fit 


Fixed n oc r 2 


Best Fit 


Fixed near 2 


Best Fit 


Fixed n oc r 2 


C f (xl(T 52 ) 


5.0 


7.7 


31 


30 


39 


32 


C T (Xl0 38 ) 


0.30 


1.9 xlO" 2 


6.2 x 10" 4 


6.8 x 10" 4 


2.4 x 10" 3 


3.4 x 10" 3 


a R 


0.9 


0.9 


0.9 


0.9 


0.8 


0.9 


a B 


-1.4 


-1.0 


-1.0 


-1.0 


-0.8 


-1.0 


s 


2.8 


2.0 


2.0 


2.0 


1.5 


2.0 


R Q (cm) 


2.6 x 10 15 


2.7 x 10 15 


5.0 x 10 15 


5.0 x 10 15 


6.7 x 10 15 


6.2 x 10 15 


vo (c) 


0.10 


0.10 


0.20 


0.19 


0.26 


0.24 


Bo (G) 


8.4 


4.2 


1.3 


1.3 


1.5 


1.6 


£ (xl0 47 erg) 


10 


2.7 


1.7 


1.8 


5.2 


5.1 


n e ,o (xlO 4 e~ 1 cm 3 ) 


170 


83 


2.2 


2.3 


1.7 


2.2 


7m,0 


1 


1 


1.8 


1.8 


3.0 


2.9 


MoCxlO^Mgyr 1 ) 


13 


6.5 


0.60 


0.63 


0.83 


0.92 


Modulation Epoch 


40,135 




1708 




474 




Abruptness Factor 


0.27,1.2 




7.7 




0.21 




Modulation Factor 


< 2.5,11 


2.9 


40 


53 


<0.20 


<0.28 



NOTE. — The subscript "0" denotes the parameter values at Af « 10 days since explosion. 



